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We report the first observation of tiie doubly Cabibbo-suppressed decays D'^ — s- K^r)^''^ using a 
791 fb"'^ data sample collected with the Belle detector at the KEKB asymmetric-energy e+e~ col- 
lider. The ratio of the branching fractions of doubly Cabibbo-suppressed relative to singly Cabibbo- 



suppressed D+ -)■ TT+ry*'^ decays are B{D+ -^ K+r))/B{D+ -> -k+t)) = (3.06 ± 0.43 ± 0.14)% and 
B{D+ -^ K+r]')/B{D+ -^ n+r]') = (3.77 ± 0.39 ± 0.10)%. From these, we find that the relative 
final-state phase difference between the tree and annihilation amplitudes in D"*" decays, 5ta, is (72 
± 9)° or (288 ±9)°. We also report the most precise measurements of CP asymmetries to date: 



= (+1.74 ± 1.13 ± 0.19)% and A^ 



+ _ 

CP 



= (-0.12 ± 1.12 ± 0.17)%. 



PACS numbers: ll.30.Hv, ll.30.Er, 13. 25. Ft, 14.40.Lb 



Decays of charmed mesons play an important role 
in understanding the sources of SU(3) flavor symme- 
try breaking structure [l|, y] and can also be sensi- 
tive probes of the violation of the combined charge- 
conjugation and parity symmetries (CP) produced by 
the irreducible complex phase in the Cabibbo-Kobayashi- 
Maskawa flavor-mixing matrix 3' in the standard model 
(SM). This SU(3) flavor symmetry structure is not well 
studied in £)+ meson decays into two-body final states 
with an 77''', since they are all Cabibbo-suppressed de- 
cays. Examples of two-body decays with an 77''' in the 
final state are the doubly Cabibbo-suppressed (DCS) de- 
cays D~^ — >■ K'^rj'^''^ and the singly Cabibbo-suppressed 



(SCS) decays D^ 



TT^rj 



(') 



The DCS decays D+ 



K^Tj^''^ have not yet been observed. The observation of 
such modes is not only intrinsically important to illu- 
minate the meson decay process but also there is gen- 
eral interest in the experimental technique of measuring 
an extremely rare decay processes with neutral particles. 
Observation of £>+ — >■ K~^r]^'^ would complete the picture 



of DCS decays for D'^ mesons decaying to pairs of light 
pscudoscalar mesons. 

In this Letter, we report the first observation of £)+ -^ 
K+Tj^'^ decays. The DCS decays D+ -)■ K+tj^'^ together 
with D+ — >■ K^TT^ can be used to measure the relative 
phase difference between the tree and annihilation am- 
plitudes {5ta), which is an important piece of informa- 
tion relevant to final-state interactions in D meson de- 
cays. Note that experimentally one is able to determine 
only the tree and annihilation amplitudes and the relative 
phase difference between them since all decays involving 
A'^ will be overwhelmed by Cabibbo-favored decays in- 
volving a K^ , with no way to distinguish between them 
because one detects only a K'g [J] . In addition, the most 
sensitive search for CP violation in £)+ — > tt^tj^'^ decays 
is reported. Observation of CP violation in Z?"*" — >■ Tr+ry^'^ 
decays with current experimental sensitivity would rep- 
resent strong evidence for processes involving physics be- 
yond the SM 0. 

The data used in this analysis were recorded at or near 



the T(45') resonance with the Belle detector [6| at the 
e"'"e~ asymmetric-energy collider KEKB [7[. The sample 
corresponds to an integrated luminosity of 791 fb~^. 

We apply the same charged track selection criteria 
that were used in Ref. [8|. Charged kaons and pions 
are identified by requiringthe ratio of particle identifi- 
cation (PID) likelihoods [8| to be greater or less than 
0.6, respectively. For kaons (pions) used in this analy- 
sis, the efficiencies and misidentification probabilities are 
approximately 87% (88%) and 9% (10%), respectively. 
For the reconstruction of the rj meson in the _D+ — > h~^ri 
decay, where h~^ refers to either 7r+ or K^, we use the 
1] — >■ 7r"''7r~7r'' mode instead of the frequently used i] — > 77 
(rijry) mode since our event selection will include stringent 
requirements on the vertex formed from charged tracks in 
the r] decay. We find that the ry — )■ 77 mode has a small 
signal to background ratio and poor 77 invariant mass res- 
olution that prohibit the final signal extraction from our 
data. To reconstruct the r]' meson in D+ — >• h~^r]' decay, 
we use the rj' — )■ Tr~^Tr~ rj^^ decay. The minimum energy of 
the 7 from the vr" or ry is chosen to be 60 MeV for the bar- 
rel and 100 MeV for the forward region of the calorime- 
ter [9| . The decay vertex of the _D+ is formed by fitting 
the three charged tracks (/i+tt+tt^) to a common ver- 
tex and requiring a confidence level (C.L.) greater than 
0.1%. For 7r° reconstruction in D~^ — > h'^i], we require 
the invariant mass of the 77 pair to be within [0.12,0.15] 
GeV/c^ and for the r] we require the invariant mass of 
the tt+tt-tt" system to be within [0.538,0.558] GeV/c^. 
In the £>+ — > /i+jy' mode, to reconstruct the daughter 
r]jj, we require the invariant mass of the 77 pair to be 
within [0.50,0.58] GeV/c^. Furthermore, in order to re- 
move a significant tt" contribution under the rj^-y signal 
peak, we reject 7 candidates as described in Ref. [lOJ. 
The TT^TT^rijj invariant mass is required to be within the 
range [0.945,0.970] GeV/c^. The momenta of photons 
from the tt" and the n^,^ combination are recalculated 
with tt" and rj mass [ll| constraints, respectively. The 
invariant mass distributions of the h^rj'^'^ system after 
the initial selection described above are shown in Fig. [T] 
where there is little indication of signal for either of the 
DCS modes. 

In order to search for Z?+ — > K^-q^'^ decays, the follow- 
ing four variables are considered. The first is the angle 
(^) between the charmed meson momentum vector, as 
reconstructed from the daughter particles, and the vec- 
tor joining its production and decay vertices [12|. The 
second variable is the isolation x^ (x?so) normalized by 
the number of degrees of freedom (d.o.f) for the hypoth- 
esis that the candidate tracks forming the charmed me- 
son arise from the primary vertex, where the primary 
vertex is the most probable point of intersection of the 
charmed meson momentum vector and the e^e~ inter- 



action region 12|. Because of the finite lifetime of D+ 
mesons their daughter tracks are not likely to be compat- 
ible with the primary vertex. The third and the fourth 
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FIG. 1: The invariant mass distributions of K^iq [K^ri') 
in the left (right) plot after the initial selection. The solid 
histograms show ■k'^t)^'' while the dashed histograms show 
K^r)^'' final states. The two inset histograms are K^rj^'' de- 
cays with enlarged vertical scales. 



variables are the momentum of the rj^^' {p„{>} ) in the lab- 
oratory system, and the momentum of the D"*" in the 
center-of-mass system (p^+). To optimize the selection, 
we maximize Csig/VA/s where esig and Mb are the signal 
efficiency and the background yield in the invariant mass 
distribution of I?+ candidates. A uniform grid of 10,000 
points in four dimensions spanned by the four kinematic 
variables described above is used to select an optimal set 
of selection requirements using Monte Carlo (MC) sim- 
ulation samples [13]. Since we use MC samples, this is 
similar to the importance-sampled grid search technique 
in Ref. IJ]. The optimal selection for the _D+ — > K^rj 
mode is found to be: ^ < 5°, xf^Q > 10, p,, > 1 GeV/c, 
and p^+ > 3 GeV/c, and for D+ -^ K+rj' is : ^ < 5° , 
y^ 

Also 



> 5, Pn' > 1.5 GeV/c, and p^+ > 3 GeV/c. The 
same selection criteria are applied to the normalization 
modes, D+ — ?► Tr+ry*^'' . Figure [5] shows the Tr+ry*^'' and 
K'^rj^'' invariant mass distributions after the final selec- 
tions used for the branching fraction measurements. Pos- 
sible structures, for example from Df — > iir+7r~7r"'"7r'' or 
£)+ -^ K^K^tt^tt'-* due to particle misidentification or 
cross-feed between 77 and 77' are studied using MC sam- 
ples; we find no indication of such background. 



A fit is then performed for D^ 



-■q 



•^'^ candidates 



and the results are shown as the top two plots in Fig. [21 
The signal probability density function (PDF) is mod- 
eled as the sum of a Gaussian and a bifurcated Gaussian 
while the combinatorial background is modeled as a lin- 
ear background. The x^/d.o.f of fits are 0.7 and 1.4, re- 
spectively. For fits to these DCS decays, we fix the width 
of the Gaussian, the two widths of the bifurcated Gaus- 
sian, and then ratio of the normalizations of the Gaussian 
and the bifurcated Gaussian to the values obtained from 
the fits to the SCS modes in order to obtain stable fits. 
The fixed widths are scaled according to the difference 
of widths observed in the signal MC samples. We ex- 
amine possible systematic uncertainties due to this later. 
The statistical significance of the signal based on the log- 



TABLE I: Yields from the data and the signal efficiencies for 
the branching fraction measurements. Errors are statistical 
only. 



Mode 



yield 



Signal Efficiency (%) 



D+ -> K+T^ 


166 ± 23 


D+ -^ K+7]' 


180 ± 19 


D+ -> TT+r; 


6476 ± 110 


D+ -> TT + V 


6023 ± 93 



1.35 ± 0.01 
1.20 ± 0.01 
1.68 ± 0.02 
1.59 ± 0.01 



likelihood ratio is 9a and more than 10a (cr represents one 
standard deviation from the background-only hypothe- 
sis) for £)+ — >• if +77 and D~^ — >• K^rj' , respectively; the 
corresponding invariant mass distributions and fits are 
shown in the lower panel of Fig. [2] The x^/d.o.f of 
fits to the K^Tj and K'^rj' final states are 0.8 and 0.9, 
respectively. In order to compute the ratio of branch- 
ing fractions of DCS modes with respect to SCS modes, 
the signal efficiencies for the selection criteria described 
above are estimated with our signal MC samples. TableU 
lists all the information used for the branching fraction 
measurements. 
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FIG. 2: The invariant mass distributions used for the branch- 
ing fraction measurements. The top two plots are for the iv'^r] 
(left) and Tv^rj' (right) final states while the bottom two plots 
are for the K'^rj (left) and K'^rj' (right) final states. Points 
with error bars and histograms correspond to the data and 
the fit, respectively. 



The dominant sources of the systematic uncertainty in 
the branching fraction measurements are the uncertain- 
ties of the parameters that are fixed in the fits to DCS 
decays, and are estimated to be 3.4% (2.1%) for the ij {ij') 
mode. These uncertainties are determined by refitting 
the data with the fit parameters varied by one standard 
deviation. Other sources include the choice of the fitting 



functions, estimated to be 2.7% (1.0%) for the rj [r]') 
mode, and the uncertainty in the PID, estimated to be 
1.1% for the both modes. A summary of the systematic 
uncertainties for the ratio of branching fraction measure- 
ments can be found in Table |TT1 The ratios of branching 
fractions are B{D+ -^ K+t^)/B{D+ -^ tt+t]) = (3.06 ± 
0.43 ± 0.14)% and B{D+ -^ K+r]')/B{D+ -^ tt+t?') = 
(3.77 ± 0.39 ± 0.10)%. We use the measurements of 
the SCS modes from Ref. [15[ to calculate the absolute 
branching fractions. Table IIIII shows the comparison of 
our branching fractions with the best present limits from 
Ref. [15| . While the measured branching fraction for the 
K~^ri mode is in agreement with the SU(3) based expec- 
tations [ij, y, the K~^r]' mode is measured to be larger, 
by approximately three standard deviations. 



TABLE II: Summary of all relative systematic uncertainties 
for the measurements of ratios of branching fractions. 



Source 




) (%) 


{B{D + ^K + r,') 
" \ S(D + ^7r + J7') 


) (%) 


PID 

Signal PDF 
Fit method 




1.1 

3.4 
2.7 




1.1 

2.1 
1.0 


Total 




4.5 




2.6 



Using the relations in Ref. [^ , which give 

|T|2 = i\A{K+i^)\'' 



\A\' 



cos 5ta 



1 



\A{K+n")\' + \A{K+T^')\ 



+ ^'M2 



- \A{K+r^)[^ 



1 



2\T\\A\ 
3 



2\A{K+r^)\' + -\A{K+^')\^ 



+ ^0\\2 



\A{K+7:^)\ 



(1) 



where T {A) is the tree (annihilation) amplitude and A 
is the specified decay amplitude, and from the recent 
branching fraction measurement of B{D~^ — ^ iir+7r°) — 
(1.72 ± 0.20) X 10-4 [H], we find that the relative final- 
state phase difference between the tree and annihilation 
in D+ decays, Sta, is (72 ± 9)° or (288 ± 9)°. 



TABLE III: Comparison of our branching fraction results to 
the present best upper limit (90% C.L.) from Ref. [I3I. The 
first and second uncertainties are statistical and systematic, 
respectively. 



Measurement 


Belle 


Ref. [15] 


B{D+ -^ K+r]) 
B{D+ -^ K+r/) 


(1.08±0.17±0.08)xl0-* 
(1.76±0.22±0.12)xl0-* 


< 1.3x10-* 

< 1.9x10-* 



For our Acp measurement in the D+ — > Tr+Ty*-'' modes, 
we re-optimize our selection by maximizing A/s'/o's where 
as is the statistical error on the signal yield Afs in the 



simulated sample. The re-optimized requirements for 
L»+ -> TT+ry decays are: ^ < 5°, xfso > 5, p^ > 1.0 GeV/c, 
and p*jj+ > 2.5 GeV/c, and for D+ -^ TT+rj' are: ^ < 5°, 
xfso > 2, Pri' > 1-0 GeV/c, andp;^+ > 2.5 GeV/c, respec- 
tively. These requirements are slightly less stringent than 
the selection criteria used for the branching fraction mea- 
surements of DCS modes. This improves the statistical 
sensitivity on Acp by around 15%. 

We determine the quantities A(jp~^^ ^ [16| by mea- 
suring the asymmetry in signal yield 



A. 



£)+^7r+?7<'' _ 



ICC 



N, 



D- 



m 



D + - 



^ A 



CP 



A 



D-- 
rcc 

FB ^ 



a: 



(2) 



where A'loc is the number of reconstructed decays. Note 
that we neglect the terms involving the product of asym- 
metries and the approximation is valid for small asym- 
metries. The measured asymmetry in Eq. ^ includes 
two contributions other than Ac p. One is the forward- 
backward asymmetry {Ap^) due to 7* — Z'^ interfer- 
ence in e^e^ — > cc and the other is the detection ef- 
ficiency asymmetry between positively and negatively 
charged pious (A'^ ). To correct for the asymmetries 
other than Acp, we use a sample of Cabibbo- favored 
D+ — )■ 07r+ decays, in which the expected CP asym- 
metry from the SM is negligible. Assuming that Aps 
is the same for all charmed mesons, the difference be- 
tween A^,^^^^"'" and AreS^'^'''' yields the CP viola- 
tion asymmetry A^p ^ ^ . We reconstruct (j) mesons 
via the K^K~ decay channel, requiring the K^K^ in- 
variant mass to be between 1.01 and 1.03 GeV/c^. This 



is the same technique as the one developed in Ref. [17[. 
In order to obtain Acp, we subtract the measured 
asymmetry for D^ — >■ (/)7r+ from that for _D+ — > Tr+ry^'' in 
three-dimensional (3D) bins, where the 3D bins are the 
transverse momentum, pJfJj, and the polar angle of the 
7r+ in the laboratory system, cos6'Jf'^, and the charmed 
meson polar angle in the center-of-mass system, cos 9* ^ . 

Simultaneous fits to the D7~, and D7 ^ invariant mass dis- 
tributions for each bin are carried out. A double Gaus- 
sian for the signal and a linear function for the back- 
ground are used as PDFs for Z?+ — > (pn^ . The average 

value over all bins is found to be A, 



DJ- 
rec 



= (0.17 ± 

0.13)%. After the subtraction of Are! ^ component, 
weighted averages of the Acp values summed over the 
3D bins are (-hi. 74 ± 1.14)% and (-0.12 ± 1.13)% for 
D"*" — ^ TT'^r] and D"*" — >■ Tr~^r]' , respectively, where the un- 
certainties originate from the finite size of the D^ — > vr+Ty 
(1.13%), D+ -^ TT+r]' (1.12%), and D+ -^ 07r+ (0.13%) 
samples. The x^/d.o.f values summed over the 3D bins 
are 28.7/11=2.6 for D+ -^ tt+t? and 15.7/11=1.4 for 
D+ -^ TT+r/. 

The dominant source of systematic uncertainty in the 



Acp measurement is the uncertainty in the A^ei de- 

termination, which originates from the following sources: 
the statistics of the Df -^ (/)7r+ sample (0.13%), pos- 
sible detection asymmetry of kaons from (p ^ K^K^ 
(0.05%) !l8| and the choice of binning for the 3D map 
(0.12%, 0.01%), for D+ ->■ TT+r] and D+ -^ tt+t;', respec- 
tively. Another source is the fitting of the invariant mass 
distribution (fit interval, choice of the fitting function), 

which contributes uncertainties of 0.05% to A^jp^^ ^ , 
and 0.07% to Afjp~^^ '' . Possible systematic uncertain- 
ties due to the fixed signal PDF parameters are estimated 
to be 0.01% for A 



D+- 



^^^ and 0.07% for A?^" 



By 
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iCP a.ii'u u.ui /o lui J-i^cP 

combining all sources in quadrature, we obtain A^jp 

= (+1.74 ± 1.13 ± 0.19)% and A^p"*''^''' = (-0.12 ± 
1.12 ± 0.17)%. These are the most precise measurements 



of A 



CP 



to date. 



In conclusion, we report the first observation of DCS 
D+ — > K+rj^'^ decays using a 791 fb^^ data sample col- 
lected with the Belle detector at the KEKB asymmetric- 
energy e~'"e^ collider. The ratios of branching frac- 
tions of DCS modes with respect to the SCS modes are 
B{D+ -^ K+T])/B{D+ -^ TT+T]) = (3.06 ± 0.43 ± 0.14)% 
and B{D+ ^ K+r]')/B{D+ ^ tt+t?') = (3.77 ± 0.39 ± 
0.10)%. Using our DCS branching fractions and that of 



D" 



K+tt" from Ref. [151, the first measurement of the 



relative phase difference between the tree and annihila- 
tion amplitudes in D+ decays is reported with 6ta = 
(72 ± 9)° or (288 ± 9)° using the technique suggested in 
Ref. [4]; this is important information relevant to final- 
state interactions. We also search for CP asymmetries 
in SCS modes down to the 0{%) level. 
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